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Abstract. 

Dark gamma-ray bursts (GRBs) are sources with a low optical-to-X-ray flux ratio. Proposed 
explanations for this darkness are: i) the GRB is at high redshift ii) dust in the GRB host 
galaxy absorbs the optical/NIR flux iii) GRBs have an intrinsically faint afterglow emis- 
sion. Within this framework, GRB 100614A and GRB 100615A are extreme. In fact, they 
are bright in the X-rays, but no optical/NIR afterglow has been detected for either source, 
despite several follow-up campaigns began early after the triggers. We analyze the X-ray 
data and collect all the optical/NIR upper limits in literature for these bursts. We then build 
optical-to-X-ray spectral energy distributions (SEDs) at the times at which the reddest upper 
limits are available, and we model our SEDs with the extinction curves of the Milky Way 
(MW), Small Magellanic Cloud (SMC), and the attenuation curve obtained for a sample of 
starburst galaxies. We find that to explain the deepest NIR upper limits assuming either a 
MW or SMC extinction law, a visual extinction of A v > 50 is required, which is extremely 
unlikely. Since both GRBs are bright in X-rays, explanation iii) also cannot explain their 
dark classification, unless optical radiation and X-rays are not part of the same synchrotron 
spectrum. An alternative, or complementary explanation of the previous possibility, involves 
greyer extinction laws. A starburst attenuation curve gives A v > 10, which is less extreme, 
despite still very high. Assuming high redshift in addition to extinction, implies an A v > 10 
at z = 2 and A v > 4 - 5 at z = 5, regardless of the adopted extinction recipe. A different, 
exotic possibility would be an extremely high redshift origin (z > 17 given the missing K 
detections). Population III stars are expected to emerge at z ~ 20 and can produce GRBs 
with energies well above those inferred for our GRBs at these redshifts. Mid- and far-IR ob- 
servations of these extreme class of GRBs can help us to differentiate between the proposed 
scenarios. 
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1. Introduction 

Long-duration gamma-ray bursts (GRBs) are 
high energy phenomena linked to the death of 
massive stars, emitting most of their radiation 
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in the hundreds of keV range. The gamma-ray 
(or prompt) event is followed by an afterglow 
at longer wavelengths, which has proven cru- 
cial to understand the physics of these sources 
and to investigate the nature of their surround- 
ing medium. 
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While the X-ray (0.1-10 keV) afterglows 
are detected in virtually all GRBs, the opti- 
cal/near infrared (NIR) counterparts are more 
elusive. Just a few months after the discov- 
ery of the first optical counterpart of a GRB 
(van Paradijs et al. 1997), the search for an 
optical afterglow associated to GRB 970828 
was unsuccessful (Groot et al. 1998) leading 
to the definition of 'dark burst' for a GRB 
with an X-ray counterpart, but not an optical 
one. Initially, this lack of detection was widely 
ascribed to the delay time between the GRB 
event and optical observations, since ground- 
based facilities could be on target only hours 
after the trigger, when the afterglow had faded 
below their sensitivity limit. 

The main scientific driver of the Swift satel- 
lite was to facilitate the GRB afterglow de- 
tection by quickly repointing its narrow field 
instruments and disseminating the GRB coor- 
dinates worldwide. The quick detection with 
XRT and the increasing number of ground- 
based automated facilities dramatically im- 
proved the optical follow-up success, reducing 
the fraction of dark GRBs. The definition of 
dark GRB was reconsidered in a more quan- 
titative manner, based on the expected spec- 
tral behaviour in the context of the most com- 
monly accepted scenario of the fireball model. 
Van der Horst et al. (2009) propose to classify 
a dark GRB using only the spectral indices of 
the power-law slopes in the optical (J3o) and 
X-ray (fix) data, where p is linked to the flux 
by the relation F v oc v~P. They noted that, re- 
gardless of many assumptions on the specific 
electron energy distribution, if both the opti- 
cal and X-ray radiation are produced by syn- 
chrotron emission from the same source, Po 
and p x are linked. In particular, Po - Px - 0.5 
if the cooling frequency lies between the op- 
tical and the X-rays, and Po = fix otherwise. 
Thus, the optical-to-X-ray spectral index al- 
lowed range is fix - 0.5 < Pox ^ Px, with 
Pox = Px-0.5 if a spectral break is present just 
below the lowest X-ray energy detected. GRB 
afterglows with p x < Px - 0.5 are classified 
as dark in this picture. 

The optical darkness can be ascribed to dif- 
ferent factors (see e.g., Perley et al. 2009): i) 
the GRB can be at high redshift, so that the 



Lyman-a absorption prevents optical identifi- 
cations; ii) dust in the GRB host galaxy or 
along the line of sight can absorb the optical 
afterglow counterpart; iii) the optical faintness 
can have an intrinsic origin. In this case, how- 
ever, also the X-ray counterpart is expected to 
be faint. 

Here we study the 'darkness' properties 
of two GRBs (namely, GRB 100614A and 
GRB 1006 15 A), which are very bright in X- 
rays, but are not detected in the optical/near-IR 
band and have no reported host galaxy candi- 
date. The next section summarizes the obser- 
vations of these two sources and our analysis 
method. Sect. 3 illustrates our results, and in 
Sect. 4 we discuss our findings and draw our 
conclusions. Decay, photon, and spectral in- 
dices are indicated with a, F, and p, following 
the standard convention t~ a , N~£, and v' 13 , re- 
spectively. 

2. Observations and analysis 

GRB 100614A and 15 A were both discovered 
by Swift-BXT and immediately repointed by 
XRT and UVOT XRT reported an extremely 
bright transient in both cases, while UVOT 
did not detect any optical counterpart. The two 
sources were followed-up by several ground- 
based facilities, beginning 4 hr and 24 min- 
utes after the trigger, respectively, but no op- 
tical/NIR counterpart down to the K band has 
been detected. GRB 1006 15 A has been ob- 
served even 9.2 days after the trigger in a deep, 
K band exposure, but again no counterpart has 
been found. All the data collected for these two 
bursts are plotted in Fig. 1. We refer the reader 
to D'Elia & Stratta (2011) for a full descrip- 
tion of the observations of these GRBs and for 
a comprehensive reference list. 

2.1. Analysis 

We extract the optical to X-ray spectral energy 
distribution (SED) for both GRBs by selecting 
those epochs at which we have the deepest and 
reddest observations (i.e. less affected by any 
dust extinction), to constrain at best the intrin- 
sic optical afterglow flux upper limit. In ad- 
dition, we attempt to select an epoch not too 
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ponent than the one responsible for the after- 
glow emission. Magnitudes reported in liter- 
ature were corrected for Galactic absorption 
(E(B-V) = 0.03 and 0.047 for GRB 100614A 
and 15A, respectively). 

A broken power-law model was fitted to 
the XRT data following the van der Horst et 
al. (2009) method, therefore fixing the SED 
normalization and the high-energy spectral in- 
dex to the value obtained from our X-ray data 
analysis (within its 90% confidence range), the 
spectral break at the X-ray energies and the op- 
tical to X-ray spectral index as fiox — fix - 0.5. 

We model the optical suppression from the 
X-ray extrapolation assuming either a Milky 
Way (MW) or a Small Magellanic Cloud 
(SMC) extinction curve. We also test the atten- 
uation curve obtained for a sample of starburst 
galaxies (Calzetti et al. 1994). We consider the 
upper limits as positive detections, and we ver- 
ify that the model-predicted fluxes are consis- 
tent with the data, i.e. equal to or below the 
upper limits. 



Table 1. Lower limits to the visual extinction 
towards GRB 100614A and GRB 1006 15 A. 
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Fig.l. The GRB 100614A (upper panel) and 
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GRB 100615A (bottom panel) X-ray lightcurves, 
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D'Elia & Stratta 2011 for details). Cyan, blue, ma- 
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which they have been obtained. Such observations 
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are used in our analysis. The faintest X-ray data in 
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close to the initial X-ray steep decay, which 
is thought to be produced by a different com- 
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Fig. 2. GRB 100614A SEDs at two different 
epochs: GTC ugriz upper limits taken 4.3 
hours after the trigger (upper panels) and the 
NIR J and K band taken 9.5 hours after the trig- 
ger (lower panels), with Pox fixed at Px - 0.5 
(solid red line) and the optical suppression 
modelled with the MW or starburst extinction 
curves. Black dotted line represents the best-fit 
value for the X-ray spectral slope (fix = 1.50). 
Dashed red lines enclose the 90% uncertainty 
in Pox- 



3. Results 

We find that GRB fluxes computed from 
both optical and NIR data are well below 
the most conservative extrapolation from X- 
rays and require strong absorption using the 
data taken within 1 day after the trigger. For 
GRB 1006 15 A, very late-time data are avail- 
able (9.2 days after the trigger), for which the 
NIR flux is still below the X-ray extrapolation. 

The GRB 100614A SED with the reddest 
flux upper limit (i.e. 570 minutes after the trig- 
ger) requires a rest-frame V-band dust extinc- 
tion of Ay > 47, 58, and 11 mag, assuming a 
MW, SMC, or a starburst attenuation curve, re- 
spectively, and using the best-fit value for fix 
(Fig. 2). Even fixing the optical to X-ray en- 
ergy spectral index to its lowest allowed value 
(within its 90% confidence range), that is, in 
the most conservative case, results still provide 



very high Ay lower limits (Table 1). Less strin- 
gent constraints on Ay are obtained using the 
GTC ugriz flux upper limits obtained 258 min- 
utes after the trigger. We obtain Ay > 13 mag 
with either the MW and SMC extinction curve 
and Ay > 8 mag with the starburst case. 

For GRB 1006 15 A, we obtain even tighter 
lower limits. The SED extracted 24 minutes 
after the burst requires Ay > 64, 79, or 16 
mag assuming either a MW or SMC extinction 
curve, or a starburst attenuation curve, respec- 
tively. These lower limits are still very high 
for the SED extracted 5.5 hours post burst: 
Ay > 58, 72, and 15 mag for the three ex- 
tinction recipes. Less critical but still high val- 
ues, are obtained even 9.2 days from the burst: 
Ay > 22, 27, and 6 mag (Fig. 3). The latter 
lower limits are lower than the ones obtained 
at earlier epochs (but still extreme), possibly 
due to a selection effect, since at later times the 
X-ray flux decreases, but the optical/NIR up- 
per limits can not become fainter consistently, 
owing to the instrument detection limits. 

Table 1 reports all the Ay lower lim- 
its evaluated from the GRB 100614A and 
GRB 1006 15 A data at the given mean observa- 
tion epochs. These lower limits are computed 
for the reported ranges of the X-ray spectral 
index (fix = F - 1, estimated in Sect. 3) that 
corresponds to the minimum, maximum, and 
mean value of the estimated 90% confidence 
range. Bold face characters indicate the results 
obtained with the most probable X-ray spec- 
tral index. Upper limits to the optical-to-X-ray 
spectral indices Pox are also shown. 

4. Discussion 

We have analyzed two Swift 'dark' GRBs, 
namely, GRB 100614A, and GRB 100615A. 
These GRBs are dark according to every def- 
inition proposed until now. They are not de- 
tected in the optical/NIR down to very faint 
limits, despite follow-up campaigns at ground- 
based facilities began minutes to hours from 
the BAT triggers (see Table 1). In addition, 
their optical-to-X-ray spectral indices satisfy 
Pox <Px~ 0.5 (van der Horst et al. 2009 crite- 
rion). The identification of these two GRBs as 
dark bursts is the consequence of their intense 
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GRB 100615A 




Fig. 3. GRB 100615A SEDs at three different 
epochs: GROND grizJHK upper limits taken 
24 minutes after the trigger (upper panels) and 
Gemini-North K and Ks upper limits taken 
5.5 hours and 9.2 days after the trigger, re- 
spectively (middle and lower panels), with fiox 
fixed at fix - 0.5 and the optical suppression 
modelled with the MW or starburst extinction 
curves. Black dotted line represents the best-fit 
value for the X-ray spectral slope (fix = 1.35). 
Dashed red lines enclose the 90% uncertainty 
infiox- 



X-ray flux coupled to the optical/NIR missing 
detections. 

The outcome of our analysis is surprising. 
To explain the deepest NIR upper limits (i.e. 
the less affected by dust) in terms of a flux sup- 
pression described by either a MW or SMC 
dust extinction laws, Ay > 47 (Ay > 58) mag 
is needed for GRB 100614A (GRB 1006 15 A) 
before one day and Ay > 22 si 9 days after 
the trigger for GRB 1006 15 A. Such extreme 



Ay values have never been observed before and 
require an explanation. 

While a SMC-like extinction curve can ad- 
equately fit a large fraction of the dust extinc- 
tion from GRB host ISM, our present picture 
of GRB host galaxies makes an extremely ob- 
scured environment of this kind a very unlikely 
possibility. Indeed, GRBs, even reddened ones, 
are hosted by blue or normal galaxies, which 
are commonly detected in ordinary galaxy sur- 
veys. This favours a scenario of a host mor- 
phology where the line of sight to the GRB 
is dusty, i.e., dust obscures only localized re- 
gions (see e.g., Perley et al. 2009 and refer- 
ences therein). An in situ obscuration appears 
to be insufficient to be responsible for the ex- 
treme extinction levels we measure. 

All these considerations hold if the optical 
radiation and X-rays are part of the same syn- 
chrotron spectrum. They could originate from 
different emission processes or even be pro- 
duced in different, independent emission re- 
gions. This is possible in particular during the 
so-called "shallow phase", where X-ray emis- 
sion may be dominated by an emission com- 
ponent that differs from the one from which the 
optical flux originates (e.g. Zhang 2007). Since 
the two SEDs of GRB 100614A, the first one 
of GRB 1006 15 A, and (marginally) the second 
one of GRB 1006 15 A are all extracted during 
the shallow phases of these GRBs, a differ- 
ent origin of the optical and X-ray emission 
in these epochs could at least in part explain 
the optical darkness of our GRBs. However, 
we note that for GRB 1006 15 A we have ex- 
tracted a SED at a very late time, about nine 
days after the end of the plateau phase, and we 
have still obtained very high Ay lower limits 
(Ay > 20) assuming either a MW or a SMC 
extinction curve. These values are less extreme 
than that obtained using the other SEDs, but 
still very high, suggesting another or at least a 
concurring mechanism to account for the GRB 
darkness. 

This complementary explanation could be 
that local extinction recipes, such as MW or 
SMC ones are inadequate for reproducing the 
optical suppression in the host galaxies of 
these two GRBs. For example, modelling the 
dust absorption using greyer extinction laws, 
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Fig. 4. The visual extinction Ay as a func- 
tion of redshift for the second epoch SED of 
GRB 1006 15 A. Green solid, red dotted, and 
blue dashed lines are for the Milky Way, Small 
Magellanic Cloud, and starburst recipes, re- 
spectively. 



such as the attenuation curve obtained from 
the observations of starburst galaxies proposed 
by Calzetti (1994), brings GRB 100614A and 
GRB 1006 15 A to require Ay lower limits that 
are less extreme, despite still being very high. 

A mixture of moderate-to-high redshift and 
extinction can reduce the dust level necessary 
to explain the SEDs. In Fig. 4, we plot as an 
example the Ay values as a function of z ob- 
tained for the second epoch of GRB 1006 15 A 
observations. Although the visual extinction is 
considerably lower than for the z = case, 
Ay ~ 10 at z = 2 and A v ~ 4 - 5 at z = 5 is still 
required, regardless of the adopted extinction 
recipe. Thus, an intrinsic origin and/or dust ex- 
tinction, coupled to a moderately high redshift 
could explain the darkness of our GRBs. 

A more exotic but intriguing possibility 
would be that these GRBs are extremely high 
redshift events. Assuming that the lack of any 
detection in the reddest NIR band (K-band) is 
due to Lya absorption from the intergalactic 
hydrogen neutral fraction, we can set a red- 
shift lower limit of z > 17. The first population 
of very massive stars (PopIII, 1O 2 M < M < 



10 3 M G ) is expected to form at z ~ 20. Then- 
death is supposed to leave behind black holes 
of several tens of solar masses, which could be 
the early progenitors of active galactic nuclei. 
These fast-spinning black holes have a rota- 
tional energy of ~ 10 55 erg or more that can 
power a GRB explosion. The isotropic ener- 
gies of GRB 100614A and GRB 100615A as- 
suming z = 18 are 1.3 x 10 54 erg and 7.2 x 
10 53 erg, respectively (using the BAT fluences). 
However, the reported isotropic energies must 
be considered as lower limits both because 
of the conservative choice of the z used and 
because the BAT detector does not constrain 
the position of the peak emission preventing a 
bolometric estimate of the emitted energy. The 
lack of detection of any host galaxy candidate 
for these GRBs represents additional support 
of this scenario. On the other hand, against the 
high redshift interpretation there is that the ex- 
pected GRB rate at z > 17 is extremely low, 
between 0.5 and 1 GRB every 10 yr (Bromm 
& Loeb 2006). 

A possible way to differentiate between ex- 
tremely high-redshift, exotic extinction recipes 
and emission from distinct components for 
sources such as these, would be to search for 
the afterglow in the mid- or far-IR bands. A 
non-detection also in these bands could hardly 
be explained using any extinction law and 
would definitely rule out a high redshift origin. 
A multiband detection compatible with the X- 
ray flux at late times (e.g. > 10 s s after the trig- 
ger, i.e. after the end of the "shallow phase") 
would instead favour the fireball model. In this 
case, given the lack of NIR detections, the ori- 
gin of the darkness would be in exotic extinc- 
tion or high redshift, depending on the spectral 
shape in the mid- and far-IR bands. 
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